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Abstract. 

We are interested in the energy-momentum relation for a moving composite in 
relativistic quantum mechanics in many-particle Dirac models. For a manifestly 
covariant model one can apply the Lorentz transform to go from the rest 
frame to a moving frame to establish an energy-momentum relation of the form 
\J (M*c 2 ) 2 + c 2 |P| 2 where M* is the kinematic mass. However, the many-particle 
Dirac model is not manifestly covariant, and some other approach is required. We 
have found a simple approach that allows for a separation of relative and center of 
mass contributions to the energy. We are able to define the associated kinematic 
energy and determine the energy-momentum relation. Our result can be expressed as 
a modified deBroglie relation of the form 



*') y/[M*(P)(?] 2 



j 

where the kinematic mass M* will depend on the total momentum P for a general 
noncovariant potential. The prefactor that occurs we associate with a time dilation 
effect, the existence of which has been discussed previously in the literature. 
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The many-particle Dirac model describes relativistic quantum mechanics of particles 
with spin 1/2. It is widely used for relativistic atomic, molecular, and nuclear structure 
calculations. For relativistic kinematics, other relativistic quantum mechanical models 
derived from quantum field theory are much more widely used since they have been 
constructed to be manifestly covariant. Although the many-particle Dirac model is 
useful for static problems, it is not used very often for kinematic calculations involving 
many-particle composites away from the rest frame. 

A potential advantage of the many-electron Dirac model is that it is much simpler 
than quantum field theory, the Bethe-Salpeter equation, and other modern relativistic 
models; we would like to be able to use it systematically to a wider class of problems 
than just static structure calculations. At issue is how the energy scales with the center 
of mass momentum, since the many-particle Dirac model is not manifestly covariant. In 
this regard, consider the discussion given recently my Marsch [lj , in which the two-body 
Dirac-Coulomb problem is solved in the rest frame; in the case of a hydrogen atom in 
motion the energy is given as 



where Eq is the rest frame energy. This is what we would expect in the case of a 
manifestly covariant theory [2]; however, the two-body Dirac-Coulomb model is not 
manifestly covariant, and it is not clear that we should expect this to come out as 
a result of a calculation done with the two-body Dirac-Coulomb model. A relativistic 
model need not exhibit a covariant energy based on the kinematic center of mass energy, 
as demonstrated by Artru [3j. This point has been emphasized recently by Javenin [I]. 

In a theory that is manifestly covariant one can take advantage of the Lorentz 
transform to establish how the energy of a composite varies with the total momentum. 
In a relativistic quantum model that is not manifestly covariant, we would not expect 
a Lorentz boost of a wavefunction in the rest frame to be consistent with a moving 
wavefunction in that model. Consequently, in general we must use some other approach 
to determine the relation between the energy for a composite at rest and in motion. 

This problem is closely related to the separation of the center of mass and relative 
degrees of freedom in relativistic quantum mechanics. We would expect that in a 
consistent scheme in which relative and center of mass degrees of freedom are separated, 
the center of mass problem should be that of a free particle in the absence of an external 




(1) 
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potential. Probably we would be concerned if this were not the case in a particular 
relativistic quantum model. 

There is a substantial literature on two-body equations (and also for more 
complicated systems) in relativistic quantum mechanics as well as quantum field theory, 
as discussed in [5], [6], [7], [H], and pi]. 

Our interest in the problem is motivated by issues rather different than those that 
have been of concern to previous authors. We have recently been interested in the 
impact of configuration mixing, or of an excited state superposition, of a composite 
on the center of mass dynamics. The basic issue here is that if the different internal 
states have different energies, then they have different masses, which should impact the 
ratio of momentum and velocity in the nonrelativistic regime. To study such problems 
we would like to use the many-particle Dirac model as a foundation if possible (since 
it is simpler than field theory, and simpler than manifestly covariant models derived 
from field theory), since we require a relativistic model to describe the different masses 
associated with different excited states. A model that is approximately correct at low 
center of mass momentum is sufficient for us, which means we do not require the model 
to be precise in the limit that the composite approaches the speed of light. Finally, we do 
not require a manifestly covariant theory, but instead we would like an understandable 
energy-momentum relation that has the correct (or at least understandable) scaling at 
low center of mass momentum. 

We found a very simple approach that allows for a separation of the center of mass 
and relative contributions to the eigenvalue in the many-particle Dirac model. The 
resulting relation between eigenvalue and energy was unexpected (at least by us). A 
part of the relation can clearly be associated with an energy-momentum relation of the 
form of Equation (CQ), which allows for a consistent definition of the kinematic mass. 
There is in addition a prefactor, which may have physical significance within the theory, 
in the context of a modified deBroglie relation. 
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Although our approach is general, it is simplest to focus on the two-body version of the 
problem here. The associated relative problem has been of interest since the early days 
of quantum mechanics [9], [ID] , and the associated time- independent relative equation 
is sometimes called the Kemmer- Fermi- Yang equation [11], |12| . 



2.1. Hamiltonian 

The two-body Dirac Hamiltonian is 

H = cti ■ cpi + a= 2 • cp 2 + /?imic 2 + f3 2 m 2 c 2 + V(r 2 - ri) (2) 
The a and (3 matrices are 

--(-) '-(::) 

where the er matrices are Pauli matrices. Only relative forces are included in this model; 
the center of mass problem here is that of a composite particle in free space. 

2.2. Classical center of mass and relative variables 

We define the classical center of mass and relative variables according to 

MR = m 1 r 1 + m 2 r 2 (4) 
r = r 2 - ri (5) 
P = P1 + P2 (6) 



P Pi _ Pi 

fi m 2 mi 

The total mass M is 



(7) 



M = m 1 + m 2 (8) 

and the relative mass \i satisfies 
111 

- = — + — (9) 

/x mi m 2 

The classical position and momentum vectors are then 

ri =R-_ r r 2 = R + -r (10) 

Pi = -P-P P2 = — P + P (11) 

m 2 mi 
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2. 3. Hamiltonian in terms of center of mass and relative variables 

The Hamiltonian can be written in terms of the quantum momentum and relative 
position operators as 

H = ai-cf-^-P-pV«2-cf— P + p)+ P 1 m 1 c 2 + I3 2 m 2 c 2 + V(r) (12) 
Vm 2 / Vmi / 

This can be recast as 



H = \M ai + ^ CK 2j-cP + (a 2 -a 1 )-cp + ^ 1 m 1 c i + /? 2 m 2 c 2 + K(r)(13) 
This is consistent with Barut and Stroebel 



2-4- Time- dependent problem and eigenvalue equation 
The time-dependent two-body Dirac equation is 

in^*(r 1>r2 ,t) = = H^(r u r 2 ,t) (14) 
We assume a stationary state solution of the form 

tt(r 1>ra ,t) = e^$( ri ,r 2 ) (15) 
where $(r l5 r 2 ) satisfies the eigenvalue equation 

M(n,r 2 ) = M(r u r 2 ) (16) 
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3. Center of mass and relative contributions to the eigenvalue 
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In the nonrelativistic problem, we are able to split the Hamiltonian cleanly into center 
of mass and relative pieces. Although there is no equivalent separation in relativistic 
quantum mechanics, it appears to be possible to develop a simple separation of the 
associated Hamiltonians (as we discuss in this section). 

We note that Barut and coworkers have discussed the separation of relative and 
center of mass degrees of freedom for this problem [13], [Hj, [15]. Of the relevant papers 
in the literature, our approach is most closely related to the ideas in these papers; 
however, our conclusions are quite different. 



3.1. Eigenvalue equation and expectation values 

Consider now the time-independent two-body Dirac equation 



mi 



m 2 



+ jjOl 2 )-cP+(ct 2 -cti)-cp+f3imic 2 +/3 2 m 2 c 2 +V(r) 



$(17) 



Suppose now that $ is an exact solution. If so, then the eigenvalue may be expressed 
in terms of expectation values according to 



Oil H Oi 2 ■ cP 

M M ' 



[a 2 - ai) ■ cp + (5imic 2 + (3 2 m 2 c 2 + V(r) 



3.2. Kinematic mass 

Suppose that we now add and subtract kinematic mass terms to the expectation values; 
this allows us to write 



Oil H ol 2 ■ cP 

M M ' 



m 2 M 



11 A + ^ ) M*c 2 



m 2 



$)(19) 



(ct 2 -oti)-cp+l3imic 2 +l3 2 m 2 c 2 +V(r)-\ 

If we assume that the contribution of the kinematic mass to the energy appears only 
through the first term, then the second term must vanish 
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[a 2 -a 1 )-cp+/3 1 m 1 c 2 +/3 2 m 2 c 2 +V(r)-( -j£0 1 + ^(5 2 ) MV 



0(20) 



This allows us to define the kinematic mass from the raio 



M*c 2 = 



<• 


(ck 2 - ai) • cp + ftmic 2 + (3 2 m 2 c 2 + V(r) 




<• 


m, m 2 
~M + ~M 


•> 



(21) 



Note that it is possible that the kinematic mass defined in this way will have different 
values for different choices of the center of mass momentum, since the two-body Dirac 
model with arbitrary V(r) is not manifestly covariant. 

3.3. Eigenvalue and energy 

With the kinematic mass defined in this way, the eigenvalue then simplifies to 



hu = ^$ 

We can rotate to obtain 
hu - (V 



> + > > + 



(22) 



(23) 

If we assume that $' has a plane-wave dependence on the center of mass coordinate 



~ e 



iP-R/ft 



then we may simplify to 

hu -. (V 



m, m 2 
~M + ~~M 



$'^(M*c 2 ) 2 + c 2 |P| 2 
We recall that the relativistic energy for a moving composite is 
E(P) = v /(M*c 2 ) 2 + c 2 |P| 2 



(24) 



(25) 



(26) 



with a fixed kinematic mass. We see that the eigenvalue hu is proportional to the 
energy-momentum relation, but a prefactor is evident. We may write 



hu(P) = ($'(P) 



mi m 2 
~M + ~~M 



$ , (P)) v /[M*(P)c 2 ] 2 + c 2 |P| 2 (27) 
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3.4- Rest frame kinematic mass 

Bound state computations are usually performed in the rest frame, which suggests that 
the kinematic mass in the rest frame (where probably the potential model will be most 
relevant) should be estimated first by solving 



and then evaluating 



(a 2 - ai) • cp + ^rnxc 2 + f3 2 m 2 c 2 + V(r) 



(28) 



* 2 



M*c 



3. 5. Extension to many particles 



<• 


(a 2 - Qi) • cp + Pimic 2 + /? 2 ra 2 c 2 + V(r) 


•> 


<• 


mi 8 + m2 8 


.) 



P=0 



(29) 



This approach can be extended to the case of many Dirac particles, as discussed in the 
Appendix. The eigenvalue in this case is of the form 



where the kinematic mass is 



£^*')V( M * c2 ) 2 + c2 l p l 2 

3 



(30) 



M*c 



* 2 



J2 a 3 ■ c *j + 2 z^-c 2 + ]r Vjh(€k - £j 

3 3 j<k 



(31) 
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From the discussion above, we have found that the many-particle Dirac equation leads 
us to a deBroglie relation of the form 



flu 



$'V(M*c 2 ) 2 + c 2 |P| 2 



In light of the energy-momentum relation from the Lorentz transform 



E(P) = A /(M*c 2 ) 2 + c 2 |P| 2 
we might view this result as consistent with a modified deBroglie relation 



&)E 



(32) 



It might be asked why such a prefactor should occur? After some thought, it seems 
that we might associate it with a time-dilation effect associated with the potential in 
the case of a bound composite. There is precedent for such an effect in the literature 
(see m, P2])- 
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This study was motivated by our interest in problems where the internal energy of 
a composite impacts the kinematics and dynamics at low center of mass momentum. 
The many-particle Dirac model is perhaps the simplest relativistic quantum mechanical 
model that might be relevant, but it has not been used much in the literature for 
problems outside of the rest frame. Since the many-particle Dirac model in general is 
not manifestly covariant, we are not able to use a Lorentz transformation to understand 
how the energy scales with momentum. 

We found a simple approach which allows us to separate the center of mass and 
relative contributions to the eigenvalue. The resulting eigenvalue relation in general can 
be written in the form 



where we have identified M* as a kinematic mass (which we would expect to be 
independent of P in a covariant many-particle Dirac model). 

If we view the eigenvalue relation as a modified deBroglie relation 



then the question to be faced concerns the interpretation of the prefactor. It seems that 
we might best associate this with a time dilation effect, consistent with recent literature 
predicting such an effect. 




(33) 
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Appendix A. Kinematic mass in the many-particle Dirac model 

The arguments given in the main text for the two-body problem can be extended to 
more Dirac particles directly, as we discuss below. 

Appendix A.l. Classical center of mass and relative coordinates 
We begin by defining the classical center of mass coordinate through 

MR = 5>,-r,- (A.l) 

j 

where 

M = (A.2) 

j 

Relative coordinates can be defined according to 

= rj-R (A.3) 

If there are N particles, then we may define N relative coordinates, but one of them is 
redundant because 

E^ = (a- 4 ) 

j 

Appendix A.2. Classical momentum coordinates 

The total classical momentum is the sum of the individual particle momenta 

P = J> (A.5) 

j 

The relative momenta are defined according to 

*3 = P, " (A.6) 
Once again there are more relative momenta than required, since 

= EPn-^E m . = ( A - 7 ) 

3 3 3 
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Appendix A. 3. Many-particle Dirac Hamiltonian 
The many-particle Dirac Hamiltonian is 

h = E«i -cp + E/W^ + E^K-^) ( A - 8 ) 

3 3 3<k 

where we presume that the potentials Vjk are functions of relative coordinates. We can 
recast this in terms of center of mass and relative coordinates to obtain 

H = E *j ■ cfo + ^PJ + E ^mjC 2 + E V jk (t k - & (A.9) 

3 3 j<k 



Appendix A. 4- Solution and expectation values 

We assume that $ satisfies the time-independent Dirac equation 

The eigenvalue then can be expressed in terms of the wavefunction as 



(A.10) 



hu — ($ 



E a i • cjx 3 + E /?i^-c 2 + e Vjkitk - ij) 

3 3 3< k 



(A.ll) 



Appendix A. 5. Kinematic mass 

We can add kinematic mass terms and write 



$)(A.12) 



E cv^+E /W+E v^*-*,-)- E mv 

Once again we require that the eigenvalue have no contribution from the relative 
problem, which occurs if we compute the kinematic mass according to 



M*c 



* „2 



E a 3 ■ C *3 + E &3 m 3^ + E V 3k(€k - £j) 

3 3 3<k 



(A.13) 
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Appendix A. 6. Eigenvalue and energy 

The eigenvalue for the many particle Dirac problem in this case is 



We can rotate to obtain 



M 



r*„2 



$' W(M*c 2 ) 2 + c 2 |Pl 2 



assuming that 
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